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It has long been recognized that active quantities 
in complex systems of many types arc not distributed 
evenly across a domain but cluster strongly into irreg- 
ular bubbles, as in a Swiss cheese. The nomenclature, 
the nature and shape of the bubbles, and the physics in 
each subject is substantially different: spottiness in high 
Reynolds number fluid turbulence HSU Si, 110, 1113 
and boundary layers |lOl | ; droplet formation in spin- 
glasses 0, 0, clustering behaviour in networks 
[liL lla : the preferential concentration of inertial parti- 
cles [la[l7l[lM[l9ll20| with applications to rain initiation 
by cloud turbulence |2ll| : clustering of luminous matter 
111 m 1^ US and magnetic bubbles in astro-physics 
|2q . are some examples. These clusters display strong 
features whose typical length scales arc much shorter 
than their averages, thus raising the question of the na- 
ture of the interface between them and the surrounding 
longer scale regions. For instance, in spin glasses the 'sur- 
face' of the droplets has a fractal-like structure whereas 
the droplets themselves have the full domain dimension 
[T3| . In fluid turbulence the concentrated sets on which 
vorticity accumulates are tubes and sheets, although the 
fractal nature of these is unclear. These sets dominate 
the associated Fourier spectra which display a spikiness 
that is the hallmark of what is usually referred to as in- 
termittency H S S, S 0, H S 0. The ubiquity of 
this irregular bubble-like topology suggests the existence 
of a set of underlying organizing principles in complex 
multi-scale systems. Using simple but broadly applicable 
mathematical ideas, this paper will demonstrate that the 
dominant physical principle behind clustering is the exis- 
tence of a hierarchy of length scales whose members are 
in competition. 

Consider a d-dimensional system whose smallest char- 
acteristic (integral) scale L is such that the system is sta- 
tistically homogeneous on boxes fl = [0,L]'^. Moreover, 




FIG. 1: An illustrative slice through Q for one value of n: 
the black kernels are surrounded by green (gray) halos. Very 
small-scale behaviour concentrates on the black & green(gray) 
regions which constitute the set At {CnKn > 1). The halos 
have scaling exponents lower than those of the black kernels. 



it is endowed with the following two properties. Firstly, 
at each point x S fi, it possesses an ordered set of length 
scales £„ = (^n{x) associated with a hierarchy of features 
labelled by n > 2 



L> £i> £2 >...>£„> i 



n+l ■ 



(1) 



The in could be thought of as an ordered set of correla- 
tion or coherence lengths; their inverses K„(a;) = £~^{x) 
clearly obey 1 < Lk„ < Lk„+i. The second assumption 
is that the ensemble averages of the iK„(x) arc bounded 
above by some ordered, positive parameters of the system 
satisfying 1 < i?„ < Rn+i 



1 < L (k„) < R„ . 



(2) 



The ensemble average (•) is a spatial average with re- 
spect to the Lebesgue measure over ft. Thus, while the 
ordering of the ln{x) must be respected at each point, 
the in themselves could be quite rough; e.g. they could 
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consist of a series of step functions. If they become very 
small near points x* then they must obey > 0{r'^~'^) 
{r — \x — x*\ and e > 0) so as not to violate 

Following an idea used in ref. [23|, consider the real 
arbitrary parameters < < 1 and < a < 1 such that 
li + a = 1. Use Holder's inequality, \AB\ < ^\A\p + ^\B\i 



where 



1, with p = n ^ and q 



< 

Re-arranging and factoring out a term (k"} gives 




Kn+l 



> «) 



(3) 



(4) 



Lower bounds on the ratio (k^) / (k„) can be found from 
(0) thereby turning into 



11"/^ 



>0. (5) 



While it is possible that the integrand in Q could 
be positive everywhere in Q, this cannot be assumed; 
the generic case is that the integrand could take cither 



sign[40|. With the definition d = L R, 
pair of inequalities 



we have the 



(6) 



for which > is valid on regions where the integrand is 
positive, designated as good regions, and negative (<) on 
bad regions. The term (£„k„)'' on the right hand side of 
© remarkably contains the arbitrary parameter p, which 
lies in the range < ^ < 1. Its existence is important 
because the ordering in ^ makes it clear from © that 
everywhere within the bad regions (<) there are large 
lower bounds on k„ with exponents containing 1/ fi 



LKyi ^ Rri 



1/p 



(7) 



Let A^i^ be the set on which > 1 and the set 

on which £„Kti < 1- Then all the bad regions (<), des- 
ignated by the clusters of black kernels in Figure 1, lie 
in A^. The green/gray halos also lie in and corre- 
spond to those parts of the good regions (>) neighbour- 
ing the bad. It is in these halos where the lower bound 
(£„K„)^ becomes operative. The white areas of Figure 
1 correspond to A~ in which the be randomly 

distributed subject to their ordering in (^). It is clear 
from (jSJ that the existence and location of the clusters 
may differ for each n. A physical picture that displays 
all clusters for every n would be the union A'^ ~ 

To show that the volume of A^ comprises a small 
part of 17, Chebychev's inequality relates the normalized 



Lebesgue measure to(^+) to the integral of iK„ over A^ 



LKn dm > m{A+) R]!^ = L"'^V+i?,V" • (8) 



Together with the relation Lk„ dm < (Lun) < Rn 
we have 



< R, 



4+1 



(9) 



Hence m{A^) is significantly smaller than unity and de- 
creases as Rn increases. Thus A^ can fill, at most, a 
small fraction of fi. With such sparse information it is 
difficult to estimate the Hausdorff' or the fractal dimen- 
sions of A^, but it is still possible to estimate scaling 
exponents |28l |. This entails making a third assumption 
of self-similarity to estimate the smallest number of balls 
of radius A+ needed to cover A^ ■ Defining A+ as 



i/p 



(10) 



for some p > 1 , it is clear that cannot be large enough 
when p = 1 because of However, any value of 
p ^ 1 will do that makes large enough to be a member 
of An- The simplest and worst estimate would be to 
write 



(11) 



Inequality (O, however, shows that A^ occupies only a 
small fraction of 51. A multiplicative factor of to(^+) is 
introduced thus 

- rn(A+,) l^-^y = {Lk+^f ■ (12) 

Instead of using ^ to estimate m{An), an assumption 
of self-similar scaling is introduced that requires that the 
change in volume of the balls with respect to n should 
scale as (the volume of A^) scales to L'^. Thus 



V+ 



K+i 



(13) 



We observe that the definition of the set A^ in principle 
involves the length scales L and A+, but not overtly X^+i- 
Yet the good and bad sets involve all three scales; L, 
A+ and A^^^. The self-similarity assumption (fT^ is an 
assumption about the nature of the set A^ that relates 
successive length scales A^ and X ^+i in an ad hoc, yet 
reasonable, fashion. Using H13|l in 112|1 we have 



A, 



n+l 



A 



Xn J (-^nfc^+i) 



2d 

T^^f/^- (14) 



From these, two estimates for JV^f emerge, one each for 
the green/gray halo and black kernel regions of Figure 1, 
whose scaling exponents |42j are independent of p 



\r+ < J {^nkn) 



(green/gray halo) 
(black kernel) 



(15) 
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TABLE I: Summary of conclusions regarding the sets A:, 
and the coloured regions in Figure 1. 



Figure 1 


black 


green (gray) 


white 


Set 








Inequality lO 


< (bad) 


> (good) 


> (good) 


Exponent 


< d 


< d{l - ft) 


= d 



For the former, the > direction of the inequality in © 
has been used together with a simple Holder inequality 

(<(i+'^))'^'>«)(^+^)/^ = (fc+)i+^ (16) 

whereas for the latter k„ < has been used. In 

contrast, without any evidence of contraction of volume, 
the formula corresponding to l(T^ for A/"" is 




where fc„ satisfies Cn^n ^ 1- The uniform scaling expo- 
nents in (|15|) are bounded by 

K^.^o < dil - m) 2?+.ke.. < d (18) 

whereas V- ^ d from The coefficients Rt^^ in lO 
to (|17|l reflect the fact that this effect is taking place only 
at length scales smaller than LRn^^^. The green/gray 
halo clearly plays the role of an interface of small but 
finite thickness between the d-dimensional (white) outer 
region and the (black) inner kernel whose dimension can 
be as high as d but could be less. When V'^^^^ saturates 
its upper bound we have 

^ihalo<rf(l-M)<^?+ker='^- (19) 

For the green/gray region to have an exponent at least 
d — 1 (a surface), ^ would lie in the range < < l/d. 
Without equations of motion, a numerical experiment 
would be necessary to estimate the i?„ by finding the 
maximum value of the ensemble average In prin- 

ciple could then be found from numerical estimates of 
^crit ^ LRn^^^ within the black kernels although if the 
K„ take very large values there it might not be possible 
to achieve resolution, /i itself may have upper and lower 
bounds that are themselves n-dependent, as in ref. p^ . 

We now proceed to discuss some examples. The first 
ideas on clustering came more than half a century ago 
from Batchelor and Townsend who observed inter- 
mittent behaviour in their high Reynolds number flow 
experiments, closely followed by observations in bound- 
ary layers by Emmons |lClj |. Batchelor and Townsend 
called this phenomenon 'spottiness' and suggested that 
the energy associated with the small scale components 
is distributed unevenly in space and roughly confined to 
regions which concomitantly become smaller with eddy 
size . Mandelbrot then suggested that these clustered 



sets on which energy dissipation is the greatest might be 
fractal in nature |30j . In measurements of the energy dis- 
sipation rate in the atmospheric surface layer, Meneveau 
and Sreenivasan interpreted the intermittent nature of 
their signals in terms of multi- fractals 01 ; a newer gener- 
ation of experiments measuring intense dissipation in tur- 
bulent fiows have been pursued by Zeff et al Q. Sreeni- 
vasan and Bershadskii have recently suggested that the 
clustering of high frequencies in a turbulent signal can be 
characterized by a scaling exponent . 

The extremely rapid time evolution of sets of high vor- 
ticity or strain in fluid turbulence is an important is- 
sue; many computations exist showing how these take 
on the nature of quasi-one-dimensional tubes and quasi- 
two-dimensional sheets which have short lifetimes [a, W, ■ 
An alternative to studying the problem in a statistical 
manner is to include time in the ensemble average (•), in 
which case the semi-infinite nature of the time-axis sug- 
gests a different measure might be necessary |4^. With 
specific reference to the Navier-Stokes equations, analysis 
is not advanced enough to deal with the full space-time 
equations (except see ref. HJ); conventional methods of 
analysis use Sobolev norms to L^-average over space and 
remove the pressure |^ [s^, 113 leaving only time. In ref. 
|27l | a hierarchy of k„ have been constructed which are 
comprised of ratios of norms (of derivatives of order n) 
and therefore functions of time only; thus the clusters 
of Figure 1 are merely gaps in the time-axis. It is then 
necessary to prove that they are finite in width and de- 
creasing with increasing Reynolds number. This involves 
finding bounds on fi. 

The second example is that of the low-temperature 
phase of spin glasses [3 113 - Our conclusions regard- 
ing the fundamental role played by the competition be- 
tween members of a hierarchy of length scales is con- 
sistent with the observation of ultrametricity in spin- 
glasses, a term that is used to denote the presence of 
a hierarchy of scales . This has been observed in 

computations on the low-temperature spin glass phases 
of the Shcrrington-Kirkpatrick |34l ISa . l37l | and Edwards- 
Anderson models , as well as in dynamic phenomena 
in complexity [1^. The results of this paper, and par- 
ticularly with reference to H19|l . are consistent with the 
droplet theory [TTl IT^ IT^ where the kernel of the droplet 
is of full dimension d but its surrounding 'surface' has a 
scaling exponent < d. In fact, Palassini and Young ^3 
have shown numerically that V^^^^ = 2.58 ± 0.02 when 
(i = 3 and Pit, i„ = 2.77 ± 0.02 when d = 4. 

In conclusion, we have shown that in a system endowed 
with a competitive hierarchy of correlation lengths, a 
clustering effect ensues in which length scales smaller 
than a critical value, and much smaller that the ensem- 
ble average scale, aggregate into small intense regions. 
The kernels of these intense regions are surrounded by 
halos that have scaling exponents smaller that of the do- 
main dimension d. We have expectations that this idea 
of competition between scales may be a useful paradigm 
in explaining the behaviour of multi-scale systems. 
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